Nonequilibrium Zeeman-splitting in quantum transport through nanoscale junctions 
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We calculate the differential conductance G(V) through a quantum dot in an applied magnetic 
field. We use a Keldysh conserving approximation for weakly correlated and the scattering-states 
numerical renormalization group for the intermediate and strongly correlated regime out of equilib- 
rium. In the weakly correlated regime, the Zeeman splitting observable in G{V) strongly depends 
on the asymmetry of the device. In contrast, in the strongly correlated regime the position Ak of 
the Zeeman-split zero-bias anomaly is almost independent of such asymmetries and of the order of 
the Zeeman energy Aq. We find a crossover from the purely spin-fluctuation driven Kondo regime 
at small magnetic fields with Ak < Aq to a regime at large fields where the contribution of charge 
fiuctuations induces larger splittings with Ak > Aq as it was observed in recent experiments. 
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Quantum transport through interacting nanoscale 
junctions has attracted much interest over the last 
decade [T]. In the Coulomb blockade regime single elec- 
tron transfer through a quantum dot coupled to two leads 
is blocked due to the charging energy U = e^/(2C) of 
the device, where C is the capacitance [5]. For quantum 
dots with an odd number of electrons, spin-flip scatter- 
ing opens up a new transport channel due to the Kondo 
effect [3]. The zero-bias conductance increases for de- 
creasing temperature and even reaches the unitary limit 
value ^ characteristic of a fully transparent device [3]. 

Recent experiments have measured the splitting of this 
zero-bias anomaly (ZBA) of the differential conductance 
G{V) in a finite magnetic fieldjH |6]. While in some 
experiments [5] the splitting was smaller than the Zee- 
man splitting in accordance with theory, in others [5] 
a crossover to larger splitting exceeding the theoretical 
predictions was observed at large magnetic field. How- 
ever, much of the theoretical calculations [3{9] have fo- 
cused on the equilibrium situation, i.e. only the linear 
response conductance. Studies out of equilibrium are 
usually either perturbative in some quantities and 
thus restricted to special parameter values (e.g. large or 
small voltages, small Coulomb interaction, etc) or em- 
ploy methods which do not properly account for the low 
temperature Fermi liquid in equilibrium^lO . Addition- 
ally, results obtained for the Kondo-model where charge 
fluctuations of the quantum dot are explicitly excluded 
might not be directly applicable to experiments, as such 
fluctuations are inevitably present in actual devices. 

In this paper we address the question how a combi- 
nation of charge fluctuations and nonequilibrium effects 
in strongly correlated systems can indeed explain the ex- 
perimentally observed discrepancies. We will pinpoint 
the crossover from spin-fluctuation to charge-fluctuation 
driven transport with increasing bias and flelds. Ad- 
ditionally, we show that the bias dependence of the 
nonequilibrium spectral function yields a splitting of the 
ZBA in the Kondo regime which is robust against asym- 



metries in contrast to a single-particle picture used in 
mean-field descriptions. 

We employ a Keldysh based conserving approxi- 
mation in the weakly correlated regime [121 [E] and 
the scattering-states numerical renormalization group 
(SNRG)[T3| in the strongly correlated regime. Both 
methods are applicable to arbitrary bias V and reduce 
to the corresponding equilibrium theory for V ^ 0. The 
SNRG being non-perturbative covers the full range of in- 
teractions but is numerically much more expensive. Its 
accuracy is determined by discretization errors [H]. Since 
the SNRG coincides with the Keldysh approach at small 
U the error remains well controlled as demonstrated in 
Ref. [13] . The temperature is set much smaller than any 
characteristic energy scale of the problem and is thus 
considered as T — >• 0. 

Since the Kondo temperature Tr- is exponentially de- 
pendent on J7[3], devices must be operated[Tl [SJ |B] in 
an intermediate regime with a moderate U to observe 
the Kondo effect experimentally. Consequently, spin and 
charge fluctuations are not well separated. The latter be- 
come significant at large bias and lead to a modification 
of the splitting as shown below. We include the charge 
fluctuations by considering the single impurity Anderson 
model (SIAM) comprising a single orbital on the quan- 
tum dot coupled to two leads, 

n= de{e- fia)cl „^c,^^^+ Y E^n^ (1) 

a,a=L,B. cr=±l 

+Un^ni + y2J— de {dlc^^^ + cl^^d^} . 

d^ and d\ (c^^^ and cl^^) are the annihilation and cre- 
ation operators for electrons on the quantum dot (in lead 
a = L, R) with spin a = ±1, ha- = d^d^, and U is the 
charging energy of the dot. We allow for different cou- 
pling Ta of lead a to the dot, but flx the total coupling 
strength, F = + Ffl, which is used as the unit of en- 
ergy. The asymmetry of the coupling is characterized by 
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the ratio R = Tji/Tl. In a finite magnetic field H the 
single-particle levels of the quantum dot are shifted by 
the Zeeman energy Aq = g^sH, E„ = Ed + f Aq. For 
simplicity, we use a constant density of states between 
—D and D for the noninteracting lead electrons in the 
wide-band hmit D U,\E^\,T a- 

The total voltage drop across the junction is given 
by the difference in the chemical potentials of the two 
leads, V = ^ fijj. In order to model the individual 
voltage drops at each contact we employ a serial resis- 
tor model, where the voltage at each contact is inversely 
proportional to its coupling strength. Keeping the dot- 
levels voltage independent, the chemical potentials in the 
leads are then given by /zl = and = —j^V. 

Two important limits are incorporated, (i) For symmet- 
ric coupling {R — 1) half the voltage drops at each side, 
and nonequihbrium effects are most pronounced, (ii) In 
the tunneling regime (i? or l/R. — )• 0) only an infinites- 
imal current flows through the junction. The quantum 
dot is in equilibrium with the stronger coupled lead and 
resembles an usual SI AM [3], where the weaker coupled 
lead acts only as a probe similar to a STM tip. 

The current / through the quantum dot is determined 
by the nonequihbrium spectral function /9o-(a;, V) 15J 

^ = ^J2j'^ [/i(^) - /«(^)] '^rp.C^, V), (2) 

where Go = x (i+ny ' ^i^^ /"(^) = i+e3(^^-M.,) are the 
Fermi functions. The differential conductance G{V) = 
dl/dV reveals the Zeeman-splitting of the quantum dot 
levels in a finite magnetic field. The peak positions in 
G{V) depend on the energetic level position of the dot, 
denoted by edot, as well as on the coupling asymmetry 
R. 

In the noninteracting case the dot spectral function 
remains voltage independent. Thus, the differential con- 
ductance at zero temperature and arbitrary R is given 
via the equilibrium spectral function [TB] 

GEq(^^) = [PEq(Mi?,) +RPEq{^J'L)] , (3) 

and maxima in G^qiV) emerge when one chemical po- 
tential coincides with a dot level, ~ Cdot- For i? > 1, 
the dominant contribution comes from the second term 
in Eq. ^ and a pronounced peak occurs for voltages 
Kiax = (1 + ;^)edot- For symmetric coupling i? = 1, 
where both terms contribute equally, a factor of two 
between peak positions in GEq{V) and PEq(w) results, 
Knax — 2edoti while in the tunneling regime R oo, 
GEqiV) (X pEq(w = V). 

This behavior extends to the weakly correlated regime. 
The differential conductance in a finite magnetic field 
Aq — 2r is shown in Fig. [T] for various values of R and 
a moderate Coulomb interaction U = —lE^ — T. For 
comparison, we added the corresponding equilibrium es- 
timation GsqiV) as thin black dashed lines. Even though 
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FIG. 1. (Color online) Normalized differential conductance 
G{V)/Go of a particle-hole symmetric quantum dot [Ed = 
— 17/2) with moderate Coulomb interaction (7 = F in a finite 
magnetic field Aq = 2r and various coupling asymmetry pa- 
rameters R as function of bias voltage V . The curves are cal- 
culated in self-consistent second Born approximation within 
the Keldysh formalism [TS]. For comparison, the crosses are 
SNRG result for 7i = 1. The estimates G^qiy) are also shown 
as thin dashed black lines for each R. The inset shows the 
position of maxima as function of R. 

there are pronounced deviations for R = \ due to the U- 
drivcn voltage dependence in V)^ there is an overall 
agreement between the equilibrium estimate GEq(V^) and 
the true nonequilibrium conductance G{V). Both results 
quickly converge upon increasing R and agree almost per- 
fectly already for i? = 10. 

Additional to the bare Zeeman splitting e^J^" = ±^ 
for a noninteracting dot, the self-energy corrections [13] 
caused by the finite U shift e^^^.^ to slightly larger val- 
ues. For the parameter values of Fig. [T] the peak po- 
sitions in pEq('^) are edot ~ ±1.25r. As visible in the 
inset, the maxima in GiV) as well as GEq(V') are lo- 
cated approximately at T^nax « (1 -I- ;^)edot, and notable 
deviations occur only near symmetric coupling. The sec- 
ond contribution from the weaker coupled lead occurs at 
14iax ~ {R+ l)edot and is visible for i? = 4 as the small 
shoulder at V ~ 6.5r. 

Fig. [2]Ja,b) shows the nonequilibrium conductance of 
a strongly correlated quantum dot in a finite magnetic 
field with particle-hole symmetric and asymmetric single- 
particle energy for various R. Only positive voltages are 
shown for the particle-hole symmetric dot of panel (a) 
due to the symmetry G{—V) = G{V). The Zeeman en- 
ergy considerably exceeds the corresponding Kondo tem- 
peratures, Tk w 0.09r for Ed = -4F and Tk w 0.2r for 
Ed = — 2r. In both cases, the tunneling regime is per- 
fectly approached upon increasing R, G{V, R — oo) cx 
Pe^'^(F) which is added to the plots as black lines. 

We clearly can distinguish two voltage regimes. At 
large voltages G{V) strongly depends on Ed as well 
as on the asymmetry R. The large-voltage maximum 
is reminiscent of the atomic charge excitation of the 
dot. Its position I^nax, as indicated by the arrows in 
Fig. ^a.,h) and plotted in Fig. ^c), is again consistent 
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FIG. 2. (Color online) Panels (a) and (b): Normalized differential conductance calculated with SNRG for large Coulomb 
interaction U = SF in a finite magnetic field, Ao = F, and for various coupling asymmetries, (a) Ed = — 4F (b) Ed ~ — 2F. 
The vertical dashed lines indicate the Zeeman-energy Aq and the different curves are offset by a constant. The small arrows 
indicate the position of maxima in G{V), which are plotted in panel (c) as function of R. A]", denotes the position at negative 
voltages of panel (b). 



with Kiax ~ (1 + ;^)edot, where edot ~ 4.2r (edot ~ 6.6r) 
for Ed = — 4r {Ed = — 2r) has been obtain from the 
equilibrium NRG spectral function p^^^{uj). 

In contrast, the position of the Zeeman-split ZBA Ak 
at low voltage (also indicated by arrows) appears to be 
almost independent of particle-hole and coupling asym- 
metry as clearly seen in Fig. [2]ja,b). Ak only slightly 
increases with increasing R [see Fig. ^c)] and is essen- 
tially given by the Zeeman-energy An « Aq. In the lin- 
ear response regime ^7it9j the spectral function is approx- 
imated to be voltage independent: p{uj,V) ~ pEq(i^)- 
This would imply that the maximum of the Zeeman-split 
ZBA in G{V) should occur at Ak ~ (l+;^)wKondo, where 
i^Kondo is the position of the Zeeman-split Kondo reso- 
nance in p-Eq{uj). As a hallmark of the Kondo effect, the 
splitting of the Kondo resonance in /9Eq(i^) is strongly en- 
hanced and approaches twice the value of a noninteract- 
ing level for large fields, WKondo ~ Ao[7H5]- The position 
of the Zeeman-split ZBA would thus evolve with R from 
Ak ~ Aq in the tunneling regime to Ak ~ 2Ao for a 
symmetric junction. Therefore, this simple picture which 
partially extents to the high-voltage behavior completely 
fails for the Kondo-correlated ZBA! 

The ZBA arises from magnetic spin- flip scattering f 17], 
where the spin of the electron on and the spin of the 
electron scattered across the dot contributing to the cur- 
rent are both flipped. In a finite external magnetic field 
such a spin- flip is associated with a finite energy cost Aq, 
which must be provided by the bias voltage for the trans- 
port channel to open, V > Aq. Hence, the asymmetry 
governs only the magnitude of the peak (and might lead 
to a merge with a charge excitation peak as in Fig. ^h) 
for negative voltages) but it does not alter its position as 
this is determined in leading order only by the difference 
in energy between the initial and final state of the dot. 

So the well-known and established picture of spin- 
flip scattering as the driving force behind the Kondo- 



effect provides a clear understanding of the asymmetry- 
independence. However, the insensitivity of the ZBA to 
asymmetries observed here as well as in experiments [SI |B] 
is highly nontrivial, if viewed in the light of the Meir- 
Wingreen formula ([2]). Obviously, a simple translation 
of maxima in the equilibrium spectra pEq^Lj), to max- 
ima in the nonequilibrium differential conductance G{V) 
fails, and a full understanding of the bias dependence of 
many-body correlations is required. The voltage depen- 
dent renormalization of the transmission is essential for a 
correct description of the conductance as already pointed 
outfS]. In the Kondo regime, our non-perturbative many- 
body calculation out of equilibrium has revealed that the 
redistribution of spectral weight in p(a;,T^)|13 with bias 
occurs in such a way that the maximum in G{V) is always 
found at Ak ~ Aq, irrespective of R and Ed- 
Fig. ^B,) summarizes the results for Ak as function of 
the applied magnetic field for small interaction U/T = 1 
as used in Fig.[T] The finite width of the ZBA determines 
the threshold field above which the splitting can be ob- 
served. Here, Kondo correlations are small, and the slope 
of Ak approaches 5(1 + at large fields. The constant 
offset indicates the Hartree shift due to a finite U. 

Fig. [3]Jb,c) displays the normalized peak position 
Ak I Aq of the ZBA for a symmetric junction with U/T = 
4, 8. In both cases, the zero-field extrapolation is consis- 
tent with the asymptotic value of the equilibrium Kondo 
regime[7], limAo-i-o ^if/^o = §■ As in equilibrium^, 
Ak/Aq increases with field, and in the strongly corre- 
lated regime [panel (c)] this increase in rather slow due 
to the presence of typical Kondo-logarithms. In contrast 
to the equilibrium Kondo regime, however, a crossover 
to a non-universal behavior Ak/ Aq > 1 occurs at fields 
of the order of the charge scale, Aq ~ ^ = 4r « 451^. 
For U ~ 4r [panel(b)], the charge fluctuations domi- 
nate the physics much earlier, and the crossover is ob- 
served already at Aq F « 3Tk- We attribute this non- 
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FIG. 3. (Color online) Position of the ZBA for particle-hole symmetric quantum dot as function of Zeeman energy for Coulomb 
interaction f/ = F (a), 4r (b), and 8F (c). Panel (a) directly displays Ak calculated with the second Born approximation for 
various R, while panels (b) and (c) show the normalized position Ak/Aq for a symmetric junction R = 1 obtained within the 
SNRG. The inaccuracies due to the numerical differentiation of I{V) are indicated by the error bars in panels (b) and (c). The 
scale at the bottom measures Ao in units of the coupling F, while the top scale is in units of the Kondo temperature Tk- 



universal Zeeman splitting at large field to the nonequi- 
librium charge fluctuations. 

We can distinguish two different regimes for strong 
correlations U > vrF: (i) At fields and voltages below 
the charge excitation scale, Kondo correlations dominate 
leading to a renormalization of Ak < Aq, (ii) at high 
fields the Zeeman split Kondo peaks merge with the Hub- 
bard side bands of p{uj, V) and the Kondo effect is de- 
stroyed. The physics is mainly driven by charge fluctu- 
ations accessible at large voltages and Ak > Ao similar 
to the weakly correlated regime. 

We have indeed reproduced the qualitative findings 
of experiments[B] and have identified the non-universal 
splitting as consequence of charge fiuctuations. How- 
ever, for C//r « 8 as estimated from experimental setups, 
the crossover in our calculations occurs at Aq/Tr- w 40 
[panel (c)] . This is an order of magnitude larger than the 
values reported in experiments where the crossover typi- 
cally occurs at Aq /Tr- « 3 — 7. We can reproduce such 
values only with a smaller Coulomb interaction U/T — A. 
This might be hint for an additional transport mecha- 
nism present in experiments, which is not accounted for 
in the simple SIAM investigated here. 

In summary, we have investigated the Zeeman splitting 
of the ZBA for weak, intermediate and strongly corre- 
lated nanodevices under true nonequilibrium conditions. 
In weakly correlated devices or generally at large bias the 
transport is driven by charge fluctuations, and the split- 
ting strongly depends on the junction asymmetry. In con- 
trast, for the Kondo-correlated ZBA present at large U/T 
our approach predicts an almost asymmetry-independent 
splitting, which is in accord with experiments and can be 
understood in a simple spin-flip scattering picture. Our 
theory explains the experimentally observed crossover [6| 
from Ak < Aq to Ak > Aq from spin to charge fluctu- 
ation driven transport, but also reveals shortcomings of 
an oversimplified description using only a single level due 
to the inaccurate scale for this crossover. 
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